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Abstract 

We consider a two dimensional magnetic Schrodinger operator on a square lattice with a spatially 
stationary random magnetic field. We prove Anderson localization near the spectral edges. We use a 
new approach to establish a Wegner estimate that does not rely on the monotonicity of the energy on 
the random parameters. 
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1 Introduction 

We consider a spinlcss quantum particle hopping on the two dimensional lattice and subject to a random 
magnetic field. This model is the magnetic analogue of the standard Anderson model with a random on-site 
potential but here the magnetic field carries the randomness in the system. The fluxes through the plaqucts 
are i.i.d. random variables. For simplicity we assume that there is no external potential. 

The main result of this paper is a Wegner estimate for the averaged density of states for the Hamiltonian 
restricted to a finite box A. More precisely, we show that the expected number of eigenvalues in a small 
spectral interval of length rj is bounded by CryjAj^. This estimate exhibits the optimal (first) power of rj, 
but its volume dependence is not optimal. Nevertheless, it can be used to prove spectral and dynamical 
localization via the standard multiscale argument. 

A similar result has been obtained earlier by Klopp et al. in [9], but under a quite restrictive condition, 
namely that in a fixed domino tiling of I? the flux on each domino is deterministically zero. This condition 
was essential for the method of [9] to work since it ensured that the magnetic field was generated by a 
stationary vector potential that could be expressed in terms of independent gauges on each domino. The 
variation of the local flux thus influenced the quantum state only on a few sites. Diagonalization of a finite 
matrix then showed that near the spectral edges the energy is a strictly monotone function of the flux. This 
monotonicity observation provided the key input for the Wegner estimate in [9] . 

If the zero flux condition is removed and we consider independent fluxes on each plaqucts, then apparently 
there is no direct monotone relation between the eigenvalues and the fluxes. Prior to our recent work [3], some 
form of monotonicity has always been used in the proofs of the Wegner type estimates in context of random 
Schrodinger operators (discrete or continuum). In case of random external potential with a definite sign, 
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the monotonicity of an eigenvalue as a function of the random couphng constants is a direct consequence of 
first order perturbation theory. For sign indefinite potentials [5] , for the random displacement model [8] and 
for random vector potentials [5, 4, 12] the monotonicity could still be extracted using some special structure 
of these models but similar ideas do not seem to apply for random magnetic fields. Note that i.i.d. (or 
stationary) random vector potentials and i.i.d. random magnetic fields represent different physical models; 
typically a stationary random magnetic field cannot be generated by a stationary random vector potential. 
We refer to Section 4 of [3] for an overview and further references. 

In a recent paper [3] we have developed a new method to prove a Wegner estimate without any mono- 
tonicity mechanism. We considered the corresponding continuous model, i.e., the Schrodinger operator in 

with a stationary random magnetic field. We proved a Wegner estimate for all energies and Anderson 
localization at the bottom of the spectrum. The magnetic field was required to have a positive lower bound 
to ensure that the current did not vanish. Furthermore, the random field had to contain modes on arbitrary 
small scales to control the large momentum regime of the current. 

In the current work we extend this approach to the discrete case, where many technical complications of 
the continuum model are absent and we can thus consider very general random magnetic fields. Apart from 
independence, the only essential condition on the fluxes is that they should be separated away from and 
TT, i.e. away from the "minimal" and "maximal" fluxes on each plaquets. This condition is analogous to the 
positive lower bound on the magnetic field in the continuous model. The proof presented here is very simple 
and it highlights the essence of our new approach introduced in [3] . 

2 Model and Statements of Results 

For any subset A C Z'^ wc introduce the Hilbcrt space ^^(A), with inner product 

The discrete magnetic Schrodinger operator will be defined on "H := £'^{7?). For the precise definitions we 
will follow the notations in [9], see also [10]. 

Let £ be the set of directed edges (arrows) in Z^, i.e., 

£ = {(x^y) ■x,y £ Z^, \x~y\ = 1}. 

For a ~ [x, y) £ £, we write a = {y, x). For a E £ we will denote by and at the first respectively second 
entry of a, i.e., a = (ai^a-t). Let F be the set of unit squares [plaquets) in 1? , i.e., 

J" = {{a;i,2:i + 1} X {a:;2,.T2 + 1} : (xi, X2) G Z^} . 

We label the elements of F by the lattice point in the lower left corner of the unit square, i.e., for x = 
{xi,X2) £ Z^ we define {xi,xi + 1} x {a;2,a;2 + 1} € F. For any fx & F we define the oriented 

boundary 

dfx = {{x,x + ei), (a; + ei, a; + ei + 62), (a; + ei + 62, x + 62), (x + 62, x)} C £, 

where ei = (1, 0), and 62 = (0, 1). 

Let T := M/ (27rZ). A function A : £ T will be called vector potential or gauge if it satisfies 

A{a) = ~A{a), Va e (2.1) 

Let 

g ■.= {AeT^ : (2.1) holds} 

be the space of vector potentials. For (x,y) S £ we wfll also use the notation A{x,y) := A{{x,y)). For a 
given vector potential A G Q, the define the differential "curl" by 

dA{f) := Ma), feF, 

a€df 
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which is a function on T with values in T. We will call the magnetic field generated by A. 

The discrete magnetic Schrodinger operator H on 1? with a vector potential A^Q y& defined by 

y&?:\x-y\ = \ 

In view of the definition of Q the operator -ff (A) is self-adjoint and its quadratic form is given by 
We have the boimd < H{A) < 8, [10], and it is well known that cr{H{0)) = [0, 8]. 

Remark. A function \ : I? ^ T defines a gauge transformation. U\^{x) :— e^^'^^^tp{x), in the sense that 
H{Ax) = U\H{A)Ul with Ax{a) := A{a) + dA(a) and dA(a) := X{ai) - \{at). It is well-known that if 
A,A£Q generate the same magnetic field, d^ = dA, then H{A) and H{A) are unitarily equivalent by 
means of a gauge transformation. Thus spectral properties of H{A) only depend on the underlying magnetic 
field. 

Next we will introduce the random magnetic field. Let = {Ljf)f^jr be a family of independent, not 
necessarily identically distributed random variables taking values in T. We assume that the distribution 
of cj/ is absolutely continuous and we denote by f/ its density function defined on T. We will denote the 
probability space by VI T-^ and the expectation w.r.t. this probability measure by E. By A^^ we shall 
denote a vector potential such that AA^ = B^^. 

The Wegner estimate will hold under the assumption that is a twice continuously differentiable function 
with uniformly bounded second derivative. Moreover, we will assume that Vf is supported away from integer 
multiples of tt. Note that not only the "minimal" flux sa is excluded, but also the "maximal" flux 
w/ « TT. To formulate this precisely, for any h € (0, 7r/2) we introduce the following subset of the torus 

Tfa :=T\{(-6,&)U(7r-6,^ + &)}. (2.2) 

Assumption A{b,D). Let Vf e C^(T) with Hw/Hc^ D and suppu/ C Tf,. 

For a rectangle A C Z^, we consider the Hamiltonian H\{A) restricted to ^^(A) as follows: for ip e ^^(A) 
and a; € A we set 

[i/A(A)V'](x) :=4V'(x)- 

\x-y\ = l 

This choice of boundary conditions is referred to as simple boundary conditions, see [7] . With respect to these 
boundary conditions we state the Wegner estimate. We remark that the proof of the localization presented 
in [7] uses this choice of boundary conditions for the Wegner estimate and for the box Hamiltonians in the 
multiscale analysis. 

We introduce the cube centered at the origin of length i e N by 

Al := {a; G : max{|2;i|, \x2\} < L}. 

We shall write = Ha^ . By xe,7j we will denote the characteristic function of the closed interval [E — 
T]/2, E + f?/2]. The Wegner estimate holds in an energy interval up to 

E„it := 4- %/8 = 1.1715... 

which is the the maximal possible value for the bottom of the spectrum, see (2.3). We now state our main 
result on the Wegner estimate. 



3 



Theorem 2.1 Suppose Assumption A{b, D) holds for some b £ (0,7r/2) and a finite constant D. Then 
for any E* < i^crit there exists a finite constant C = C(b, D, E*) such that for any E > and r] > with 
E + r]/2 <E*, we have 

ETrx£;,^(HL(A„)) <C77L«. 



Remark. By using the unitary transformation tpi^) ~^ (— have the unitary equivalence 
H{A) = 8 — H{A). Therefore the Wegner estimate holds likewise for E—i^/2 > 8 — E* and the same symmetry 
applies for the result about localization (Theorem 2.2 below). The restriction E* < i^crit originates from the 
proof of Lemma 4.2. Since the Lifshitz asymptotics (2.4) only holds for small energies, a Wegner estimate 
for energies below -Ecrit is sufficient for the Anderson localization near the band edge which is our main 
application of Theorem 2.1. Nevertheless it is an interesting question on its own to establish the largest 
possible upper threshold E* of the range of energies E for which one could extend Lemma 4.2 and hence the 
Wegner estimate. 

We now state the result on localization. For this part, we assume that the random fluxes w/ are not 
only independent but also identically distributed with density v. Under this condition, 7J(A^) is an ergodic 
operator, and its spectrum a{H{A^)) = S is actually independent of w a.s., see [1] (Proposition V.2.4.). 

For E gM., the integrated density of states is defined by 

k{E) ^ lim — 5— # {eigenvalues of i^A I, (^w) < -E), 

where the limit exists u a.s. and is independent of the choice of sample w, see Appendix C of [10]. The density 
of states is independent of the choice of boundary conditions, which can be seen using the min-max principle 
[10]. Lifshitz asymptotics is shown in [9] under the assumptions that suppw C T \ (— c, c), ±c G suppw for 
some < c < TT and v is Lipshitz continuous on T \ (— c, c). Under these assumptions 

[i;o,8-£;o], with £;o = -Bo(c) :-4(l-cos(c/4)), (2.3) 

[10, 9]. In [9] it is shown (Theorem 1.2) that 

log(-log(fc(£;)) 

lim sup p F ^ < -1- 2.4 

EiEo yog{E - Eo) 

This result implies, roughly speaking, 

< as EiEo. 

for any 5 > 0. Lifshitz tail estimate and Theorem 2.1 imply localization at the bottom of the spectrum using 
standard arguments, see [2, 11, 7] for details. Here we only state the result noting that the assumptions 
ensure that Eq < Scrit and thus a Wegner estimate always holds at the bottom of the spectrum. 

Theorem 2.2 Suppose Assumption A(b, D) holds for some b G (0, 7r/2) and a finite constant D. Assume 
that the random fluxes ujf are independent, identically distributed and ±b G suppv. Then Anderson local- 
ization holds near the bottom of the spectrum. Namely, there exists an E\oc > EQ{b) such that H{Ai^) has 
dense pure point spectrum on [£'o(^'), -Bloc] almost surely, and each eigenfunction associated to an energy in 
this interval decays exponentially. 



3 Proof of the Wegner estimate: Theorem 2.1 

Let a/, denote a vector potential of a magnetic field of flux 1 through the square /. I.e., da/ = <5/ where 5f 
is the function on T which is one at / and zero otherwise. Later in Section 4 we will choose specific gauges, 
see (4.2) below. In the current section we only need that the absolute value of a/ is bounded by one. 
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Let A = Al for some i e N. Define J"a := {/ € J" : / C A} and := {a € £^ : a € A x A}. Henceforth 
we restrict the functions af to £\. We define the vector potential 

then dAuj — . Moreover we will drop the subscript lo in the notation and set H = Hl {A) in this section. 

We will use perturbation theory, wc abbreviate Yf = , and for any state if) we introduce the expectation 
of the l/-derivativc of the energy in state if): 

\x-y\ = l 

Let A be a non-degenerate eigenvalue of H with a normalized eigenfunction ip. In that case the eigenvalue 
A is a function of the random variables {w/}- For each / G J^, we have by the Hellmann-Feynman theorem 
and a straight forward calculation that 

>/A=^ = (>/i?)^. (3.1) 

To give a physical meaning to the right hand side we introduce the current of a wavefunction in the 
Hilbert space. For cp <S ^^(A), we define the current J^p as the function on £\ given by 



J^ia) := -2Re(/j(ai)ie'^("V(at)- 

The current does not depend on the chosen gauge and observe that 

J^{a) = ^J^{a). (3.2) 

If if is an eigenvector then it is a straightforward calculation to show that the "divergence" of J^p vanishes, 
namely for all x G A 

J^(a;,x + e) = 0. (3.3) 

e: |e| = l 

We can write the right hand side of (3.1) as 

{yfH)i, ^\T. «/(«)^^(«)- (3-4) 

aS£A 

Notice, since we sum over directed edges we have a factor i. We also remark that the left hand side is gauge 
invariant, while the right hand side seems to depend on the gauge a/. The gauge independence of the right 
hand side follows from the fact that the divergence of is zero (3.3). 

In Section 4 we will prove the following key technical estimate: 
Proposition 3.1 Let b G (0,7r/2) and E* < -Ecrit- Then there exists a constant C — C{h,E*) such that 

J2 {YfH)l > C-^L-^ (3.5) 
for any normalized eigenunction ip with eigenvalue in (— oo,i?*] and for any collection of fluxes ujf G T^. 
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Proof of Theorem 2.1. Let Ai, A2, ... denote the eigenvalues of H. First we assume that they are simple. 
We fix E* < Ecrit- In the sequel we set x = XE,rj, with E + r]/2 < E*. From Proposition 3.1 and (3.1) it 
easily follows that 

Trx(H) = ^x(A.) <CL^Y.T. i^fXifxiX,), (3.6) 
To estimate the right hand side wc introduce the following functions, 

Fix) := / xWdt, Giy) /' Fix)dx. 



For any Y — Yf and any £ wc have by the chain rule and Leibniz 

(YXefxiX,) = Y^G{Xe) - {Y^Xe)F{X,). (3.7) 

To estimate the sum of the second term over all eigenvalues Xi, we will use the following Lemma, for a proof 
see [3]. 

Lemma 3.2 We have for any Y = Yf 

Tr {Y^H)F{H) < Y,{Y^X,)F{Xe). (3.8) 
e 

Thus combining (3.6), (3.7) and (3.8), we have 

Tr x{H) < CL'^ (Tt YfG{H) - Tr {YfH)F{H)^ . (3.9) 

To estimate the second term we use that ||i^||oo < Grj. Moreover, we use that ||y^?i/|| < 4, which can be 
seen using that for any Lp G ^^(A) we have 

{YfHv){x)^- a}{x,y)e^^^^'y^^{y). 

\x-y\ = l 

To estimate the the first term in (3.9) we integrate by parts after taking expectation, 

ETr y/G(iJ)= / Yl z;c(wc)d^C^TrG(i?)= / [] wa^c)d^C / «}'(c^/)dc^/Tr G(i?). 

Using that G{y) < Crjy we can estimate |TrG(i7)| < CL^rj. Thus wc obtain 

ETrx(i^) < G?7i^ (3.10) 
and the Theorem follows in the non-degenerate case. 

In the argument above we used that A = is a simple eigenvalue since we implicitly assumed that Af 
is difFerentiable w.r.t. w/. Differentiability of the eigenvalues may be ensured even in case of degeneracy by 
choosing an appropriate branch, but this choice may depend on / hence a more careful treatment is needed 
which we explain now. 

Suppose that A is a degenerate eigenvalue with degeneracy m\ for some fixed lo. Let P\ denote the 
eigcnprojcction and wc fix {^p\h\h=i,...,mx orthonormal basis of the corresponding cigcnspace. Moreover, 
for each fixed / G J^a we can also choose orthonormal cigenfunctions '4'f,\,hj with /i = 1, ...,mx, which are 
real analytic functions of (if the other variables of tu are kept fixed) and their eigenvalues Ef x ^ are real 
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analytic in as well, see [6] (Theorem 1.10. Section II). By analytic perturbation theory, the derivatives 
YfEf^xjt (/i = 1, 2, . . . , m\) are the eigenvalues of the matrix 

T -.^ PE{YfH)PE \raiiPE. 

By Jensen's inequality, we have 

rnx mx mx mx 

Y^iYfEfA^)' = TrT^ = Y.i'PxM^T'vx^) > J2ivx,h,T^XM)' = ^(F/i?)^, (3.11) 

h=l h=l h=l h=l 

Summing (3.11) over /, we find, using Proposition 3.1, that for A < E*, 

mx 

^ J2iYjEf,,j,)' > mxC-'L-\ (3.12) 
/ej^A h=i 

We can now proceed as before. Using (3.12) and (3.7), we find 
Trx(i/)= "^^^(^) 

mx 

< E E T.(YfEf,xj.fx{Ef,x,h) 

Xe<y(H) feJ^A h=l 

mx 

= E EE (^/ G(i?AA,.) - iYj'Ef,,,,)FiEj,,,,)) 
feJ^A xecr{H) h=i 

<CL^ (ttY^G{H) ^TT(Y^H)F{H)y 

where we used Lemma 3.2 in the last step. The rest of the proof is the same as in the case of non-degenerate 
eigenvalues and we thus proved Theorem 2.1. |— I 

We remark that the possible degeneracy of the eigenvalues may also be treated by standard perturbation 
theory for almost all lj, using the fact that the interior of the set of w's where there is no eigenvalue crossing 
has full measure. 

4 Proof of Proposition 3.1 

We recall that by (3.4) we have 

The following lemma inverts this linear relation and expresses the current in terms of (YfH)^. 
Lemma 4.1 For a Cz £, let fa be the unique square in T such that a G df. Then for a G £a, we have 

Ma) = CaiYf^H)^ - CaiYf^H)^, (4.1) 
where Ca — 1 if fa ^ J' a and Ca — otherwise. 
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Proof of Lemma 4.1. We will introduce four different gauges, a*^"^-*, t = 1, 2, 3, 4, for the magnetic field 5f 
whose flux is 1 through the square / and zero elsewhere. For f = fx (z T with x = (.ti, X2) we set 

a, (y, y + 02) = , a, (y, y + ei) = I ^ ^^j^^^^.^^ , (4.2) 

2)^,„^„N_n ..(2). r 1 , if 2/i>xi,2/2 = a;2 



(y,J/ + ei) =0 , a)'{y,y + e2) 



, otherwise 



a/ l2/,?/ + e2j-U, (2/,y + eij - I Q ^ otherwise 

(4)/ , N n (4)^ , ^ / ^1 ' if yi < 1/2 = 2:2 

(y, y + ei) = , a, (y, y + 62) = | ^ ^ ^^j^^^^.^^ 

and extend the definition to £ by (2.1). 

Now wc can prove Lemma 4.1. Let a = (y, y + ei) G be a horizontal edge. For —L < y2, the gauge 
T = 1 gives (4.1). For the boundary case y2 = —L wc can use the gauge r = 3. The identity for the edge a 
follows using (3.2). The vertical edges follow similarly using the gauge t — 2,4:. q 

The next lemma gives a lower bound on the current of an eigenfunction. The proof will be given in 
Section 5. 

Lemma 4.2 Let b e (0,7r/2) and E* < i^crit- Then there exists a constant C = C{b,E*) such that 

aeSA 

for all normalized eigenf unctions ip of H\{A) , with energy E < E* , and for any vector potential A with 
dA e Tb (see {2.2)). 

Note that the flux is required to be separated away from and tt, i.e. not only the zero flux is excluded 
but also the "maximal" flux. The first exclusion is obvious since if dA = 0, then the eigenfunction can 
be chosen real and then clearly = 0. If dA = tt on each plaquct, then one can choose a gauge A with 
A{x,x + ei) = and A{x,x + e2) — ttxi and extend it by (2.1). In that case the Hamiltonian H\{A) is 
real. In particular, the eigcnfunctions are real for the maximal flux, hence the current vanishes in this case 
as well. 

Proof of Proposition 3.1. Using Lemma 4.1 we have 

I J^(a)P < 2ca{Yf^H)l + 2c^{Yf^H)l 

Observing that each square gives rise to 8 directed edges, after summing up this inequality for all a £\, 
we find 

32 J2 iYfH)l > E 

The proposition now follows from Lemma 4.2. |— I 

5 Proof of the regularity lemma: Lemma 4.2 

Let ^p E £'^{A) be an eigenfunction of H\{A) with energy E. If ip{x) ^ 0, we write 

^(x) = e*^(^)|V(x)| (5.1) 
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for some real function A(x) S T. Then the current is 

J^(a) ^ 2|^(a,)||^/'(at)l siii(^a), (5.2) 

with ^(^j) + A(af) — \(ai). The goal will be to find a unit square Q £ J-\ such that infa-^g is 

bounded from below. Using that the magnetic flux through Q is separated away from zero and tt, one can 
then show that there is a nonvanishing current along the boundary of that square. 
To this end let xq G A be a point where "0 attains its maximum absolute value, 

|^(a;o)| =M:=sup,gAl^(a;)|. (5.3) 
In this section we use the convention that ip{x) = if .t ^ A. 

Lemma 5.1 Let the energy E of the eigenfunction tjj satisfy < E < 4. Then the following statements 
hold. 

(a) We have 

max \^P{y)\> {1~ E/4)M. 

\y~Xa\=l 

(h) Suppose yo A is a nearest neighbor of xq with |'0(j/o)| < eM for some e > 0. Then 

min \^j;iy)\ > {2 - E - e)M. 

J/eA\{yo} 
\y-xo\ = l 

(c) Suppose 2/0 G A is a nearest neighbor of Xq with |'0(j/o)| ^ '^A/ for some k > 0. Then 

{(4 - E)k -2}M< |^(yo + t)| + myo - t)|, 

where t is a unit vector orthogonal to yo — xq. 

Proof. Adjusting the phase of the eigenvector we can assume that ipixo) = M. Since the statements of 
the lemma are gauge invariant, for notational simplicity we may choose for (a) and (b) a gauge, such that 
A{xo,y) = for all y with |xo — y\ = 1. The statement (a) follows from the eigenvalue equation 

J2 ^{y)^{i~E)iJj{xo). 

1/eA 

\v-xo\ = l 

and by taking real parts 

4 max Re'ipiy) > (4- E)'ip(xo)- 

\y-xo\ = l 

For part (b), we fix j/ G A \ {ya} such that \y — xq\ = 1. Then by the eigenvalue equation 

V'(y) = (2 - ii;)^(xo) - ^(yo) + 2V(xo) - ^ ^(z). (5.4) 

2eA\{yo/y} 

xo-z| = l 

Now taking the real part of both sides we find using that Re['0(a;o) — '(/'(•z)] > 0, 

RcV'(y) > (2 - E)^{xq) ~ eM. 

This yields (b). To prove (c), we choose a gauge such that A(j/o, = for all |yo ^ y| = 1. By the eigenvalue 
equation at yg 

(4 - £;)V'(2/o) = V'(a;o) + ^'(Syo - a;o) + '/'(yo + 1) + ^(yo - t). 
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Hence by the triangle inequality 

(4 - £;)|VXyo)| < 2Af + \^{y^ - t)| + |V;(yo + t)|, 
and (c) now follows. q 

With the above lemma we can show the following proposition. 

Proposition 5.2 Let E* < i?crit = 4 — -v/S- There exists a positive number c = c{E*) such that for any 
eigenfunction with eigenvalue E < E* there exists a cube Q £ J- a such that 

min^eQ|'i/;(a;)| > cM. (5.5) 

Proof. Let e = 1/10 and recall the definition of M and a;o from (5.3). We consider the following two cases. 

Case 1: minj,gA:|y-xo|=il'0(y)l > <^M. 

By Lemma 5.1 (a) there exists a nearest neighbor j/q G A of xq such that 

myo)\ > il-E*/4)M >0. 
Then by Lemma 5.1 (c) (and using the notation introduced there) there exists a cr £ { — 1,1} such that 

1^(2/0 + ^t)| > i {(4 - E*)il - E*/4) - 2} M. 

Thus choosing the square Q = {a;o,a;o + (7t,yQ,yQ + crt}, the estimate (5.5) follows from the choice of 
E* <4-V8. 

Case 2: mmyizA:\y-xo\=iH{y)\ < <^M. 

Let j/o e A be the nearest neighbor of xq with |?/;(2;o)| = miny£A:|j;-a;o|=ilV'(y)l- Let yo = 2xq — yo- Then 
by Lemma 5.1 (b) and (c) there exists a cr e { — 1, 1} such that 

myo+<7t)\ >^{{i-E*){2-E* -e)-2}M. 

Thus choosing the square Q ~ {xq, xq + crt, y^, y^ + crt}, the estimate (5.5) follows from Lemma 5.1 (b) and 
the choice of E* and e. q 

With the above Proposition we can show Lemma 4.2 using that a nonzero flux produces a current. 

Proof of Lemma 4.2. Since ip is £^-normalized, clearly, |'(/'(a;o)| > 1/i, where xo is defined by (5.3). By 
Proposition 5.2 there exists a unit square Q £ Ta such that for a constant c (depending on i?*), we have 

inf \i^{x)i^{y)\ > cL-^. 

x,yeQ 

With respect to the parameterization (5.1) we have by (5.2) 

^ \Ma)\'>c'L-^ ^ sm'i^a). 

aedQ aedQ 

This will imply the lemma using that for some positive constant c (depending on b) 

sm\^a) > c. (5.6) 

a£dQ 

To show this, we first note that luq — X^aeag ^ '^b- We will show that this implies that there is an 
a £ dQ such that ipa £ Tb/s- The estimate (5.6) will then follow since sin^(-) is bounded from below on T^/g 
by a positive constant (depending on b). Suppose that there would not be an a G dQ such that ipa £ Tb/s. 
Then ipa = Uai^ + ba with G Z and |6o| < b/S. In that case ujq = nir + b, with n £TL and |&| < 6/2. This 
implies loq ^ T;,, which is a contradiction. |— I 
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